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Abstract 



The ALICE experiment has measured low-mass dimuon production in pp collisions at s/s = 7 TeV 
in the dimuon rapidity region 2.5 < y < 4. The observed dimuon mass spectrum is described as a 
superposition of resonance decays (r),p,co,rj ,0) into muons and semi-leptonic decays of charmed 
mesons. The measured production cross sections for CO and (j) are < p t < 5 GeV/c,2.5 < 

y < 4) = 5.28±0.54(stat)±0.50(syst)mb and < p, < 5 GeV/c,2.5 < y < 4) = 0.940 ± 
0.084(stat) ±0.078(syst) mb. The differential cross sections d 2 a/dydp t are extracted as a function 
of pt for co and (j). The ratio between the p and CO cross section is obtained. Results for the cj> are 
compared with other measurements at the same energy and with predictions by models. 



*See Appendix [Alfor the list of collaboration members 
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1 Introduction 

The measurement of light vector meson production (p,oo,<p) in pp collisions provides insight into soft 
Quantum Chromodynamics (QCD) processes in the LHC energy range. Calculations in this regime are 
based on QCD inspired phenomenological models [CO that must be tuned to the data, in particular for 
hadrons that contain the u, d, s quarks. The evolution of particle production as a function of ^/s is 
difficult to estabilish. Measurements at mid-rapidity in pp collisions at the beam injection energy of 
the LHC (y/s = 0.9 TeV) were performed by the ALICE experiment [;2], and compared with several 
PYTHIA |j3] tunes and PHOJET HJ. The comparison showed that, for transverse momenta larger than 
~ 1 GeV/c, the strange particle spectra are strongly understimated by the models, by a factor of 2 for 
and 3 for hyperons, with a smaller discrepancy for the (p. Extending the measurements to larger energies 
and complementary rapidity domains is needed in order to further constrain the models. 

Moreover, light vector meson production provides a reference for high-energy heavy-ion collisions. In 
fact, key information on the hot and dense state of strongly interacting matter produced in these collisions 
can be extracted measuring light meson production. 

The ALICE experiment at the LHC can access vector mesons produced in the rapidity range 2.5 < y < 4 
through their decays into muon pairs Q. In this Letter we report results obtained in pp collisions at 
= 7 TeV in the dimuon transverse momentum range 1 < p t < 5 GeV/c based on the full data sample 
collected in 2010 with a minimum bias muon trigger. The measurement is done via a combined fit of the 
dimuon invariant mass spectrum after combinatorial background subtraction. 

2 Experimental setup 

The ALICE detector is fully described elsewhere |5j. The main detectors relevant for this analysis are 
the forward muon spectrometer, which covers the pseudo-rapidity region — 4 < r\ < —2.5, the VZERO 
detector and the Silicon Pixel Detector (SPD) of the Inner Tracking System. 

The elements of the muon spectrometer are a front hadron absorber, followed by a set of tracking stations, 
a dipole magnet, an iron wall acting as muon filter and a trigger system. 

The front hadron absorber is made of carbon, concrete and steel and is placed at a distance of 0.9 m from 
the nominal interaction point (IP). Its total length of material corresponds to ten hadronic interaction 
lengths. The dipole magnet is 5 m long and provides a magnetic field of up to 0.7 T in the vertical 
direction which gives a field integral of 3 Tm. 

The muon tracking is provided by a set of five tracking stations, each one composed of two cathode 
pad chambers. The stations are located between 5.2 and 14.4 m from the IP, the first two upstream of 
the dipole magnet, the third in the middle of the dipole magnet gap and the last two downstream. The 
intrinsic spatial resolution of the tracking chambers is ~ 100 /im in the bending direction. 

A 1.2 m thick iron wall, corresponding to 7.2 hadronic interaction lengths, is placed between the tracking 
and trigger systems and absorbs the residual secondary hadrons emerging from the front absorber. The 
front absorber together with the muon filter stops muons with momentum lower than 4 GeV/c. The 
muon trigger system consists of two detector stations, placed at 16.1 and 17.1 m from the IP. Each one is 
composed of two planes of resistive plate chambers (RPC), with a time resolution of about 2 ns. 

The SPD consists of two cylindrical layers of silicon pixel detectors, positioned at a radius of 3.9 and 
7.6 cm from the beam. The pseudo-rapidity range covered by the inner and the outer layer is |t] | < 2.0 
and \r] \ < 1.6, respectively. Besides contributing to the primary vertex determination, it is used for the 

'in the ALICE coordinates, the muon spectrometer covers the pseudo-rapidity range —4 < 7] < —2.5, where the z axis is 
oriented along the beam direction. However, since in pp collisions results are symmetric with respect to y = 0, we prefer to 
drop the negative sign when quoting the rapidity values. 
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input of the level-0 trigger (LO). 

The VZERO detector consists of two arrays of plastic scintillators placed at 3.4 m and -0.9 m from the 
IP and covering the pseudo-rapidity regions 2.8 < r\ < 5.1 and —3.7 < r\ < —1.7, respectively. This 
detector provides timing information for the LO trigger and has a time resolution better than 1 ns, thus 
giving the possibility to reject beam-halo and beam-gas interactions in the off-line analysis. 



3 Data selection and analysis 

During the pp run in 2010, the instantaneous luminosity delivered by the LHC to ALICE ranged from 
0.6 x 10 29 to 1.2 x 10 30 cm s . The fraction of events with multiple interactions in a single bunch 
crossing was less than 5%. The data sample used in this analysis was collected using the muon trigger, 
which is activated when at least three of the four RPC planes in the two muon trigger stations give a 
signal compatible with a track in the muon trigger system. To evaluate the integrated luminosity (JS^ n t)> 
a minimum bias (MB) trigger, independent of the muon trigger, was collected in parallel. It is activated 
when at least one out of the 1200 SPD readout chips detects a hit or when at least one of the two VZERO 
scintillator arrays has fired, in coincidence with the arrival of bunches from both sides. 

The integrated luminosity was determined by measuring the MB cross section Omb and counting the 
number of MB events. The Gmb value is 62.3 mb, and is affected by a 4% systematic uncertainty. It 
was obtained measuring the cross section CVoand 10, for the occurrence of coincident signals in the two 
VZERO detectors (V0AND) in a van der Meer scan [7]. The factor OVoand/<7mb was obtained as the 
fraction of MB events where the L0 trigger input corresponding to the V0AND condition has fired. Its 
value is 0.87 and is stable within 0.5% over the analyzed data. The full data sample used for this analysis, 
amounting to an integrated luminosity of approximately 85 nb _1 , was used to extract the p t distributions. 
Part of the data was not collected with the MB trigger in parallel with the muon trigger. For this fraction, 
the integrated luminosity could not be measured and the ft) and cross sections were determined with 
the remaining subsample corresponding to Jz^ nt = 55.7 nb _1 . 

Track reconstruction in the muon spectrometer is based on a Kalman filter algorithm (U |9]. Straight 
line segments are formed from the clusters on the two planes of each of the most downstream tracking 
stations (4 and 5), since these are less affected by the background coming from soft particles that emerge 
from the front absorber. Track properties are first estimated assuming that tracks originate from the IP 
and are bent in a uniform magnetic field in the dipole. Afterwards, track candidates starting in station 

4 are extrapolated to station 5, or vice versa, and paired with at least one cluster on the basis of a % 2 
cut. Parameters are then recalculated using the Kalman filter. The same procedure is applied to the 
upstream stations, rejecting track candidates that cannot be matched to a cluster in the acceptance of 
the spectrometer. Finally, fake tracks that share the same cluster with other tracks are removed and 
a correction for energy loss and multiple Coulomb scattering in the absorber is applied by using the 
Branson correction (SJ. The relative momentum resolution of the reconstructed tracks ranges from 1% 
at 20 GeV/c to 4% at 100 GeV/c. 

Muons were selected requiring that the direction and position of each muon track reconstructed in the 
tracking chambers match the ones of the corresponding track in the trigger stations. A cut on the muon 
rapidity 2.5 < < 4 was applied in order to remove the tracks close to the acceptance borders. Muon 
pairs were selected requiring that both muons satisfy these cuts. Approximately 291,000 opposite-sign 
{N + -) and 197,000 like-sign (N++, N ) muon pairs passed these selections. 

The opposite-sign pairs are composed of correlated and uncorrelated pairs. The former constitute the 
signal, while the latter, coming mainly from decays of pions and kaons into muons, form the combi- 
natorial background, which was evaluated using an event mixing technique. The distribution obtained 
was normalized to 2R^JN ++ N , where N ++ (N ) is the number of like-sign positive (negative) pairs 
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integrated in the full mass range. It is assumed that the like-sign pairs are uncorrelated. The fraction 
of correlated like-sign pairs, coming from the decay chain of beauty mesons and B — B oscillations iflOl 
was determined from the measured open charm content and the ratio between open beauty and charm 
(see below). It amounts to s» 0.5% for p t > 1 GeV/c and M < 1.5 GeV/c 2 , and was thus neglected. The 

R factor is defined as yJA ++ A , where A ++ (A ) is the acceptance for a ++ ( ) pair, and 

takes into account possible correlations introduced by the detector. It was evaluated using two meth- 
ods. The first employs MC simulations to determine the acceptances A±±. The other method uses the 



mixed-event pairs to estimate R as R = /^mked^y'^mixed^mixed^ w h ere ^ymked is the num b e r of mixed 
pairs for a given charge combination. The two methods are in agreement for p t > 1 GeV/c. We obtain 
R = 0.95 for pi > 1 GeV/c. The event mixing procedure was cross-checked by comparing the results 
obtained for like-sign mixed pairs with the non-mixed ones. The shapes are identical, while the number 
of like-sign pairs estimated with the event mixing is lower than the one in the data by 5%. We take this 
value as the systematic uncertainty on the background normalization. The signal-to-background ratio 
for p t > 1 GeV/c is about 1 at the (j> and ft) masses. Alternatively, the combinatorial background can be 
evaluated using only the like-sign pairs in the non-mixed data, and calculating for each AM mass bin the 

quantity 2R(AM)y^N ++ (AM)N (AM). Figure \T\ shows the invariant mass spectrum for opposite sign 

muon pairs in different p t ranges, together with the combinatorial background estimated with the event 
mixing technique or using the like-sign pairs. It is seen that the two techniques are in good agreement 
for p t > 1 GeV/c. For lower pair transverse momenta both methods fail in describing the background. 
In this region, the method based on the like-sign pairs gives a backgound mass spectrum that overshoots 
the opposite-sign pair spectrum, while the event mixing technique does not reproduce the non-mixed 
like-sign pairs spectra. The analysis is thus limited to p t > 1 GeV/c. The event mixing technique is 
used, since it is less affected by statistical fluctuations. 

After subtracting the combinatorial background from the opposite-sign mass spectrum, we obtain the raw 
signal mass spectrum shown in Fig. |2] The mass resolution at the mass is Om ~ 60 MeV/c 2 , in good 
agreement with the Monte Carlo simulation. The processes contributing to the dimuon mass spectrum 
are the light meson (tj, p, ft), rj , 0) decays into muons and the correlated semi-leptonic open charm and 
beauty decays. The light meson contributions were obtained performing a simulation based on a hadronic 
cocktail generator. The input rapidity distributions for all particles are based on a parametrization of 
PYTHIA 6.4 results obtained with the Perugia-0 tune [11]. The same procedure is followed for the 
v\' p t distribution, while for p, ft) and (j) the transverse momentum is described with a power-law function, 
used also by the HERA-B experiment to fit the pi spectrum lfl2l : 



The parameters n and po were tuned iteratively to the results of this analysis. The p t distribution of tj 
is based on preliminary results from rj production yields measured in the two-photon decay channel by 
ALICE |[L3l . The open charm and beauty generation is based on a parameterization of PYTHIA (H. 
The detector response for all these processes is obtained with a simulation that uses the GEANT3 lPT4l 
transport code. The simulation results are then subjected to the same reconstruction and selection chain 
as the real data. The invariant mass spectrum is fitted with a superposition of the aforementioned 
contributions. The free parameters of the fit are the normalizations of the TJ — )• jU/xy, ft)— > MJ., <p — > [ipL 
and open charm signals. The processes T] — > /J./J. and ft) — > are fixed according to the relative 

branching ratios. The contribution from p — > jU/i was fixed by the assumption that the production cross 
section of p and ft) are equal ITT31 IT6l IT711 . The tj' contribution was set fixing the ratio between the r\' 
and T] cross sections according to PYTHIA. The ratio between the open beauty and open charm was 
fixed according to the results from the LHCb Collaboration |[T8l [191 . The main sources of systematic 
uncertainty are the background normalization and the relative normalization of the sources, mainly due 
to the error on the branching ratios for the ft) and 7]' Dalitz decays. The raw numbers of and p + ft) 
resonances obtained from the fit are M aw = (3.20±0.15) x 10 3 and N7£ m = (6.83 ±0.15) x 10 3 . 
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Figure 1: (Color online) Invariant mass spectra for opposite-sign muon pairs in pp at y/s = 7 TeV in different p t 
ranges. The combinatorial background, evaluated from opposite-sign pairs in mixed events (red line) or like-sign 
pairs in non-mixed events (blue points), is also shown. 

4 Results 

The <p production cross section was evaluated in the range 2.5 < y < 4, 1 < p t < 5 GeV/c through the 
formula: 



N: 



MB 



9 A^BR^^l+nNMBN^ 8 ' 

where TV™* is the measured number of mesons, A§ and are the geometrical acceptance and the 
efficiency respectively, Nmb is the number of minimum bias collisions, Omb is the ALICE minimum bias 
cross section in pp collisions at ^/s = 7 TeV, and Njf B /N^ MB is the ratio between the number of single 
muons collected with the minimum bias trigger and with the muon trigger in the region 2.5 < < 4, p t > 
1 GeV/c. The number of minimum bias collisions was corrected, as a function of time, by the probability 
to have multiple interactions in a single bunch crossing. Finally, BR((j) — > = (2.95 ±0.03) x 10~ 4 is 
the branching ratio into lepton pairs. Assuming lepton universality, this number is obtained as a weighted 
mean of the measured branching ratio in pL + pL~ with that into e + e~, because the latter has a much smaller 
experimental uncertainty than the former [20]. The number of mesons was evaluated by performing a 
fit to the mass spectrum for each Ap t = 0.5 GeV/c interval in the transverse momentum range covered 
by the analysis. The acceptance-corrected results were then summed in order to obtain the total number 
of mesons. In this way the dependence of the acceptance collection on the input p t distribution used 
for the Monte Carlo simulation becomes insignificant. Alternatively, a fit was performed on the mass 
spectrum integrated over 1 < p t < 5 GeV/c and a global correction factor was applied. The results of 
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Figure 2: (Color online) Dimuon invariant mass spectrum in pp at s/s = 7 TeV after combinatorial background 
subtraction for p x > 1 GeV/c (triangles). Light blue band: systematic uncertainty from background subtraction. 
Red band: sum of all simulated contributions. The width of the red band represents the uncertainty on the relative 
normalization of the sources. 

the two approaches agree within 3%. The first approach was used for the results reported in this paper. 
The meson acceptance and efficiency correction in the range covered by this analysis was evaluated 
through Monte Carlo simulations and ranges from 10% to 13%, depending on the data-taking period. 
The ratio N^ B MB strongly depends on the data taking conditions and was evaluated as a function 
of time. 

We obtain Cty(l < p t < 5 GeV/c,2.5 < y < 4) = 0.940 ± 0.084(stat) ±0.078 (syst) mb. The systematic 
uncertainty results from the uncertainty on the background subtraction (2%), the branching ratio into 
dileptons (1%), the muon trigger and tracking efficiency (4% and 3% respectively), the minimum bias 
cross section (4%) and the ratio N^ B /N^~ MB (3%). The first two contributions have been described 
above. The others are common to all analyses in the dimuon channel, and are extensively discussed 
elsewhere |2T1 . Here, only the main points are briefly summarized. The muon trigger efficiency was 
estimated measuring the number of J/y mesons decaying into muons, after efficiency and acceptance 
corrections, in two ways: in the first case both muons were required to match the trigger, while in the 
second only one muon needed to fulfill this condition. The tracking efficiency was evaluated starting 
from the determination of the efficiency for individual chambers, computed by taking advantage from 
the redundancy of the tracking information in each station. The same procedure was applied to the data 
and to the Monte Carlo simulations. The differences in the results give the systematic uncertainty on 
the tracking efficiency. The error on the minimum bias cross section is mainly due to the uncertainties 
in the beam intensities [22] and in the analysis procedure adopted for the determination of the beam 
luminosity via the van der Meer scan. The error on the ratio N^ B /N^ was evaluated comparing the 
value measured as described above with the information obtained from the trigger scalers, taking into 
account the dead time of the triggers |[23l . 

Table Q] compares the present measurement with some commonly used tunes of PYTHIA [3 ] (Perugia- 
OH, Perugia- 1 1 G4), ATLAS-CSC GJ] and D6T E6J) and PHOJET 0. It can be seen that Perugia-0 
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Figure 3: Top: Inclusive differential production cross section d 2 o^ /dydp t for 2.5 < y < 4. The error bars repre- 
sent the quadratic sum of the statistical and systematic uncertainties, the red boxes the point-to-point uncorrelated 
systematic uncertainty, the blue box on the left the error on normalization. Data are fitted with Eq. (G3 (solid line) 
and compared with the Perugia-0, Perugia- 1 1, ATLAS-CSC and D6T PYTHIA tunes and with PHOJET. Bottom: 
Ratio between data and models. 

and Perugia- 1 1 underestimate the (j) cross section (by about a factor of 2 and 1.5, respectively), while the 
others agree with the measurement within its error. 

The differential cross section d 2 c§ / dydp t is shown in Fig. [3] (top). Numerical values are reported in 
Table |2] /? t -dependent contributions to the systematic uncertainties, due to the uncertainty on trigger 
and tracking efficiency and background subtraction, are indicated as red boxes. The uncertainty on the 
minimum bias cross section, branching ratio and Nff B /N$ MB ratio contribute to the uncertainty in the 
overall normalization. As stated above, the <p cross section is extracted from a subsample of the data 
used to determine the p t distribution, and is thus affected by a larger statistical uncertainty, resulting in a 
5% contribution to the normalization error. Fitting the expression in Eq. £T|) (solid line) to the differential 
cross section gives po = 1.16 ±0.23 GeV/c and n = 2.7 ±0.2. The PYTHIA and PHOJET predictions 
are also displayed in Fig. [3l where the bottom panel shows the ratio between the measurement and the 
model predictions. PYTHIA with the ATLAS-CSC and D6T tunes reproduces the measured differential 
cross section, while the others predict a harder p t spectrum. 

The results are compared to measurements of — > K + K~ for 2.44 < y < 4.06 by the LHCb Collabora- 
tion E71 in Fig. [4] The observed shapes of the p t distributions are similar. In order to compare with our 
integrated cross section result, the differential cross section measurement by LHCb was integrated for 
Pt > 1 GeV/c and scaled by a small correction factor, obtained from PYTHIA (Perugia-0), to account 
for the slight difference in rapidity acceptance. The result is c§ = 1.07 ±0.15(stat. ±syst.) mb. When 
the statistical errors and the part of the systematic uncertainty which is not correlated among the two 
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Figure 4: (Color online) Top: Inclusive differential production cross section d 2 o^/dydp t , as measured via the 
decay into dimuons (black triangles). The blue box on the left represents the error on normalization. The data 
are compared to the measurements in the kaon decay channel by LHCb (black open circles) f27]|- Bottom: Fit to 
the differential cross section measured in dimuons divided by the cross section measured in the kaon channel by 
LHCb. 



experiments are properly taken into account, the two measurements are in agreement. 

The ratio N$/(N P +N m )=BR(<j) ->■ iin)a$/[BR{p ->■ jUjU)a p + BR(co ->■ fifi)a a ], corrected for accep- 
tance and efficiency, was calculated for 1 < p t < 5 GeV/c, giving 0.416 ± 0.032(stat.) ± 0.004(syst.). 
Systematic uncertainties are due to the normalizations of ft) — > /I/I7T , v\' — > and combinatorial 
background. The corresponding ratio is calculated with PYTHIA and PHOJET. All the predictions un- 
derestimate the measured ratio, as reported in Table [T] The p t dependence of this ratio is shown in Fig. [5] 
The Perugia-0, Perugia- 11 and D6T tunes systematically underestimate this ratio, while PHOJET cor- 
rectly reproduces the data for p t > 3 GeV/c, and ATLAS-CSC is in agreement with the measurement for 
Pt > 1.5 GeV/c. 

In order to extract the ft) cross section, the p and ft) contributions must be disentangled, leaving the p 
normalization as an additional free parameter in the fit to the dimuon mass spectrum. The result of the fit 
for 1 < p t < 5 GeV/c gives o p /a m = 1. 15 ± 0.20(stat) ± 0. 12(syst) , in agreement with model predictions, 
as shown in Tabled] The systematic uncertainty was evaluated changing the normalizations of the r\' — > 
HHjand ft) — > pLpLTt according to the uncertainties in their branching ratios and the background level by 
±10%, which corresponds to twice the uncertainty in the normalization. The ft) production cross section, 
calculated from this ratio, is a ffl (l < p t < 5 GeV/c,2.5 < y < 4) = 5.28 ± 0.54(stat) ±0.50(syst) mb. 
This value is in agreement with the Perugia-0 PYTHIA tune, while the other tunes and PHOJET 
overestimate the ft) cross section, as shown in Table [TJ 

In Fig. [6] (top) the ft) differential cross section is shown. Numerical values are reported in Table |2] A fit 



10 



The ALICE Collaboration 



3 1 



0.8 



0.6 
0.4 
0.2 



— i — | — i — i — i — i — i — i — i — i — i — | — i — i — i — i — | — i — i — r 

Data ALICE pp \(s = 7 TeV 

PYTHIA Perugia-0 2.5 < y < 4 

PYTHIA Perugia-11 

PYTHIAATLAS-CSC 

PYTHIA D6T 

PHOJET 



p t (GeV/c) 



Figure 5: Ratio N$ / (N p +N m ) as a function of the dimuon transverse momentum. 



Table 1: Measured cross sections and ratios compared to the calculation from PYTHIA with several tunes and 
PHOJET in the range 1 < p t < 5 GeV/c, 2.5 < y < 4 



Qfr (mb) CJ m (mb) J^+nZ Qp/ff m 

ALICE llll measurement 0.940 ±0.084 ±0.078 5.28±0.54±0.50 0.416 ±0.032 ±0.004 1.15 ±0.20±0.12 

PYTHIA/Perugia-0 0.50 5.60 0.22 1.03 

PYTHIA/Perugia-11 0.62 7.81 0.20 1.03 

PYTHIA/ATLAS-CSC 0.91 6.50 0.35 1.05 

PYTHIA/D6T 1.12 9.15 0.30 1.04 

PHOJET 0.87 6.89 0.30 1.08 



of Eq. (Q} (solid line) to the data gives po = 1.44 ±0.09 GeV/c and n = 3.2 ±0.1. As shown in the same 
figure (bottom), all the PYTHIA tunes reproduce the p t slope, while PHOJET gives a slightly harder 
spectrum. 



5 Conclusions 

Vector meson production in pp collisions at yfs = 7 TeV was measured through the dimuon decay channel 
in 2.5 < y < 4 and p t > 1 GeV/c. The inclusive (j) production cross section 0^(1 < p t < 5 GeV/c,2.5 < 
y < 4) = 0.940 ± 0.084(stat) ±0.078(syst) mb was measured with a sample corresponding to an inte- 
grated luminosity J^j nt = 55.7 nb _1 . Calculations based on PHOJET and PYTHIA with the ATLAS-CSC 
and D6T tunes give results that are in agreement with the measurement, while the Perugia-0 and Perugia- 
1 1 PYTHIA tunes underestimate the cross section by about a factor of 2 and 1.5, respectively. The ratio 
N^/iNp +N m ), calculated for 1 < p t < 5 GeV/c, gives 0.416 ±0.032 ±0.004. This value is reproduced 
by PHOJET for p t > 3 GeV/c, and by the ATLAS-CSC tune for p t > 1.5 GeV/c, while the other tunes 
underestimate the ratio in the full range 1 < p t < 5 GeV/c. By measuring the ratio of the p and ft) cross 
sections, o p /o m = 1.15 ±0.20(stat) ±0.12(syst), it was possible to extract the inclusive ft) production 
cross section a ffl (l < p { < 5 GeV/c,2.5 <y < 4) = 5.28 ±0.54(stat) ±0.50(syst) mb. While all models 
correctly reproduce the measured o p / o a ratio, the ft) cross section is correctly reproduced only by the 
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Figure 6: (Color online) Top: Inclusive differential (0 production cross section d 2 a a /dydp t for 2.5 < y < 4. The 
error bars represent the quadratic sum of the statistical and systematic uncertainties, the red boxes the point-to-point 
uncorrected systematic uncertainty, the blue box on the left the error on normalization. Data are fitted with Eq. (Q]) 
(solid line) and compared with the Perugia-0, Perugia-11, ATLAS-CSC and D6T PYTHIA tunes and PHOJET. 
Bottom: Ratio between data and models. 



Table 2: <j) and CO differential cross sections for 2.5 < y < 4. Statistical, bin-to-bin uncorrelated and correlated 
systematic errors are reported. 

p t (GeV/c) d 1 ^ /dydp t (mb/(GeV/c)) d^g/dydpt (mb/(GeV/c)) 

[1,1.5] 0.695 ±0.079 ±0.046 ±0.053 3.69 ±0.35 ±0.24 ±0.32 

[1.5,2] 0.268 ±0.032 ±0.018 ±0.020 1.75±0.15±0.12±0.15 

[2,2.5] 0.147±0.014±0.010±0.011 0.857 ±0.069 ±0.057 ±0.075 

[2.5,3] 0.0665 ±0.0074 ± 0.0044 ±0.0051 0.339 ±0.029 ±0.022 ±0.030 

[3,3.5] 0.0403 ±0.0044 ± 0.0027 ±0.0031 0.220±0.019±0.011 ±0.019 

[3.5,4] 0.0169±0.0031 ±0.0011 ±0.0013 0.0880 ±0.0088 ±0.0058 ±0.0077 

[4,4.5] 0.0131 ± 0.0022 ± 0.0009 ± 0.0010 0.0648 ± 0.0062 ± 0.0043 ± 0.0056 

[4.5,5] 0.0069 ±0.0017 ±0.0005 ±0.0005 0.0301 ±0.0039 ±0.0020 ±0.0026 



Perugia-0 calculation, and overestimated by the others. The differential production cross sections of CO 
and (j> were measured. The p t dependence of the cross section agrees well with other measurements 
done in the kaon decay channel. The ATLAS-CSC and D6T tunes correctly reproduce the p t spectrum, 
while the other calculations predict harder spectra. PHOJET predicts also a slightly harder p t spectrum 
for the CO, while PYTHIA provides slopes which are closer to the one obtained with this measurement. 
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